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We study the effect of a random Flory-Huggins parameter in a symmetric diblock copolymer melt which is 
expected to occur in a copolymer where one block is near its structural glass transition. In the clean limit the 
microphase segregation between the two blocks causes a weak, fluctuation induced first order transition to a 
lamellar state. Using a renormalization group approach combined with the replica trick to treat the quenched 
disorder, we show that beyond a critical disorder strength, that depends on the length of the polymer chain, 
the character of the transition is changed. The system becomes dominated by strong randomness and a glassy 
rather than an ordered lamellar state occurs. A renormalization of the effective disorder distribution leads to 
nonlocal disorder correlations that reflect strong compositional fluctuation on the scale of the radius of gyration 
of the polymer chains. The reason for this behavior is shown to be the chain length dependent role of critical 
fluctuations, which are less important for shorter chains and become increasingly more relevant as the polymer 
length increases and the clean first order transition becomes weaker. 



INTRODUCTION 

Copolymer systems, i.e. macromolecules built of se- 
quences of chemically distinct repeat units, are of particular 
interest due to the phenomenon of microphase separation and 
the resulting formation of complex ordered structures chang- 
ing their macroscopic, mechanical properties fH S I3. 

7]. For example, a di -block copolymer melt consisting of 
blocks of A and B monomers is chiefly characterized by the 
Flory-Huggins parameter 

X = ^ab-^(vaa + vbb) (1) 

which characterizes the segregation strength between A and B 
monomers. Here, v„./ is a measure for the short distance re- 
pulsion between s and i' monomers = A or B) in units 
of k^T. For large enough % phase separation into A-rich 
and B-rich regions occurs|2|. These regions are limited in 
size due to the covalent bond between the blocks resulting 
in microphase separation. In case of equal length of the 
A and B blocks the transition between the disordered, i.e. 
mixed phase and a microphase separated state is via a fluc- 
tuation induced first order transition|4], originally proposed 
by Brazovskii|8, 9]. At low temperature the system devel- 
ops lamellar order with a period determined by the radius 
of gyration of the polymer The observation of this transi- 
tion in poly(ethylene-propylene)-poly(ethylethylene) di -block 
copolymers of approximately equal block volume is proba- 
bly the most convincing experimental verification of the Bra- 
zovskii scenario of weak crystallizationi 10|. 

The role of disorder in copolymers has been the subject 
of a number of theoretical studies. In many cases the disor- 
der is caused by fluctuations in the sequence length of the 
blocks O O Il3 H Q- It was shown that fluctuation 
effects beyond mean field theory are crucial to stabilize mi- 
crophase separation for disordered sequences 1 16, 17, 18|. In 
the limit where the typical segment length is small compared 
to the total length of the polymer, recent numerical simula- 



tions found strong deviations from theories based upon an ef- 
fective Ginzburg-Landau theory, while the latter is applicable 
if the block length, even though random, remains large with 
high probability! 19|. Using a replica theory to analyze the 
quenched random segments, only replica symmetric solutions 
were found indicating that disorder does not result in subtle 
aging and non-ergodicity behavior as found for example in 
spin glasses. A much stronger impact of quenched random- 
ness is caused by "random field" disorder This can be re- 
alized in case of a copolymer melt in a gel matrix ||2Q] with 
preferential adsorption of one of the monomers on the gel and 
was studied by Stepanow et al. fl?] who found a glassy state 
for intermediate values of % while a lamellar state only occurs 
for larger %. The dramatic impact of random field disorder be- 
comes also evident from the results of Refs.|22,'23','27] where 
it was shown for a model in the same universality class that in- 
finitesimal random field disorder leads to a state with one step 
replica symmetry breaking, i.e. a self generated glass with 
non-ergodic aging behavior 

Another realization of quenched disorder, intermediate be- 
tween random sequences and random field disorder of a gel 
matrix, occurs in case of fluctuating interaction strength be- 
tween the monomers, i.e. Vjj/ — > Vjj/ (x) which leads to a ran- 
dom, spatially varying Flory-Huggins parameter 

X^X(x) =5C + 5x(x)- (2) 

In Ref. iz^ it was argued that this is the case in a copolymer 
where one block is near its structural glass transition. The 
decreased mobility of the monomers of this block implies a 
partial annealing and leads to slow non-thermal compositional 
fluctuations which can be described by Eq|2l Systems where 
one of the two blocks are close to a glass transition were stud- 
ied for example in Refs.|26, 27]. Using a variational approach 
an ordered and a glass like state with strong compositional 
fluctuations were identified, depending on the strength of the 
fluctuations 1 25 1 . 

In this paper we analyze symmetric block copolymers with 
random Flory-Huggins parameter, Eq|2| using a renormaliza- 
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tion group approach. In agreement with Ref.f25i we find 
that a critical strength of the fluctuations 5x (x) is needed to 
destroy the fluctuation induced first order Brazovskii transi- 
tion to a lameflar state. In addition, in the glassy state we 
find, at low energies, a strong renormalization of the distri- 
bution function, P(5x(x)), characterized by strong nonlocal 
fluctuations and determine a phase diagram for variable poly- 
mer length, A^, mean Flory-Huggins parameter % and fluc- 
tuation strength, g, of 5x(x). We demonstrate that for suf- 
ficiently large chains a glassy behavior becomes inevitable, 
caused by the very different chain length dependence of ther- 
mal and disorder related fluctuations. The renormalization 
group technique used in this paper is a generalization of the 
approach of Hohenberg and Swift|28] to the case of systems 
with quenched randomness. 

Our results demonstrate how disorder affects a first order 
transition and are thus of relevance beyond the specifics of 
random copolymer systems. The fluctuation induced first or- 
der transition studied in this paper is weak. An excess entropy 
is generated by a large phase space of fluctuations, render- 
ing the system unable to reach a critical point while stabiliz- 
ing an ordered structure at a finite but rather large correla- 
tion length|8, 29]. In the opposite case of a strong first or- 
der transitions the influence of quenched disorder of the type 
discussed here was studied in Ref.|30]. In this and subse- 
quent studies! 3ji l32ll it was shown that for d = 2 disorder 
leads to a vanishing latent heat, a result obtained via mapping 
the two coexisting phases with sharp interfaces onto a ran- 
dom field Ising model. For d >2 the first order transition sur- 
vives below a certain disorder strength. The extension of the 
above mapping to d > 200 

demonstrates that the behav- 
ior for larger disorder is again in the same universality class 
as the random field Ising model. The situation is different in 
the present model. Due to the large interface width between 
distinct ordered regions at the weak first order transition, the 
mapping onto the random field Ising model cannot be justified 
anymore. Even though it is not possible for us to specify the 
precise nature of the corresponding transition in our case, our 
results strongly suggest that the transition between a regime 
with effectively weak disorder to a disorder dominated behav- 
ior is not of the random field Ising type. On the other hand, 
the fact that the resulting first order transition is weak allows 
us to obtain a qualitative understanding of our results using 
a modified Harris criterion|35|, typically used to analyze the 
role of quenched disorder at continuous phase transitions. 

The remainder of the paper is organized as follows: In the 
next section we introduce the model and summarize the main 
steps of the replica trick to perform the disorder average. We 
introduce the renormalization group approach in section III, 
where we also summarize the main results of this calculation, 
while the details of the derivation and solution of the flow 
equations are given in appendices. The results of our calcula- 
tion and a generalized Harris criterion are given in section IV 
and we briefly summarize our findings and their implications 
in section V. 



MODEL AND REPLICA TRICK 

We consider Np polymer chains (« — 1, ...,Np), each con- 
sisting of statistical segments (i = 1, ...,N). The A block 
on each molecule contains Nf statistical segments and the B 
block has N{1 — f) segments. We restrict the discussion to 
the case / = 1/2. The relevant degrees of freedom of the 
polymer are the positions of the segments R,,,;. The system is 
characterized by a Gaussian statistical weight and additional 
excluded volume pseudo-potential: 

+ ^ L v,y(R„,.)5(R„,.-R„.y). (3) 

where a is the characteristic persistence length of the polymer 
and po = is the monomer density. In the spirit of Flory- 
Huggins lattice theory the monomer density for a dense poly- 
mer melt is roughly given by po — a^^ . v^j is dimensionless 
and characterizes the excluded volume interaction. We con- 
sider a symmetric block copolymer with 
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which are assumed to vary in space. The mean repulsion be- 
tween different blocks, vab, is taken to be lai-ger than vaa and 
vbb such that for the mean Flory-Huggins parameter holds 
X > 0, i.e. A and B monomers segregate on the average. 
Following the procedure first outlined by Leibler|2] (for a 
recent review see Ref.^^), one introduces a monomer num- 
ber density fields PaW = L^j<'v 5(r — R„j) and pe (r) = 

§(r — R„,j). Using the incompressibility assumption 



of a dense melt. 



Pa (r) + PB (r) =po, 



one obtains a theory for a collective field, (|) (r), that describes 
the microscopic fluctuations of A and B monomers : 



^(r) 



c / Pa (r) 1 
y/a\ Po 2 

^ f pB (r) 1 

Vfl V Po 2 



with c=1.1019|2]. The effective action of the problem is the 
Brazovskii model|81] for a scalar order parameter (|) but with 
random mass term: 



(4) 
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Following Refs.0 [loll . To, Fq and qo can be expressed in 
terms of the parameters of the polymer Hamiltonian, Eq|3l 



to 
qo 

To 



■ iXs - X) 



1.945 

156.56 

c'^aN ' 



(5) 
(6) 

(7) 



with = 10. 49 /A'. Rg = ay ^ is the radius of gyration of a 
Gaussian chain. A natural dimensionless bare constant of the 
problem is Fo/^o = \J putting the system in the weak cou- 
pling regime for > 10^, similar to the conclusion obtained 
within in the self-consistent Hartree approachP,li,fT^. 
For the probability distribution of t (x) = -^^5% (x) we as- 



sume a Gaussian form 



;')((x'; 



(8) 



with f(x)f(x') = 2A(x — x'). In what follows we will assume 
uncorrelated disorder with correlation function: 

A(x-x') = Ao5(x-x') . 

Aq can be expressed in terms of the mean square fluctuations 

of the Flory-Huggins parameter, Ax^, as 

23.14Ax2 



Ao = 



aVN It 



(9) 



' /v fi^'x 5x(x)5x(0) characterizes the spatial 



Here, Ax^ = jy^ 

fluctuations of the Flory-Huggins parameter within the vol 



ume of a single chain Vb — coiTesponding dimen- 

sionless strength of the disorder is Aq /<7o = 4.85 ^ and does, 

Xs 

for fixed not decrease with A^. The important ratio 



Ao/Fo = 0.218^A?i/2^ 



(10) 



therefore grows for larger chain length, A^. The relative 
strength of the disorder becomes larger for longer chains. 
We determine the averaged free energy via the replica trick 
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F = -riogZ = T lim - Z" 

m^O m 



1) 



and obtain 

= |DfP[f]|D'>exp|^-p£5[^„]j 

= yD"'(j)exp(-P5W[(^]), 
with replicated Hamiltonian: 

m A f 

a=l ^ a,b-' 

In the remainder of this paper we will analyze this replicated 
action using a renormalization group approach. 



(11) 





2Aa/b 




Figure 1: Momentum shell integration procedure. The configurations 
with higher wavelength modulations are integrated out, reducing the 
low energy degrees of freedom to configurations closer to the perfect 
lamellar arrangement (as illustrated in the lower panel) 



RENORMALIZATION GROUP APPROACH 

The crucial difference between the Brazovskii model, Eq|3 
and an ordinary (|)^-model of Ising-type ferromagnets is the 
nonlocal term ( [V^ + ^q] (|)) , which strongly prefers mo- 
menta q with |q| = as opposed to the state q = 0. It is 
this enhancement in the phase space of low energy fluctua- 
tions that causes the mentioned fluctuation induced first or- 
der transition [gj. In case of ordinary n-vector (|)'*-theories 
with quenched disorder a renormalization group approach was 
developed in Refs.fsi, "3?, "3?, "3?]. However the dramatic 
change in the low energy phase space of the Brazovskii model 
requires a different formulation. 

The low energy modes of the problem are momentum states 
on a sphere with radius and one must adapt the usual dec- 
imation of high energy states in the momentum shell renor- 
malization group accordingly. Instead of the conventional 
shell integration l40bl , where momenta with Ae^' < |q| < A 
are eliminated, one has to eliminate states in a shell with 
distance A from the sphere with radius qo. This procedure 
requires A ^ ^o- This decimation approach is sketched in 
Fig^and was used in Ref.|28] to study the clean Brazovskii 
model. Here we generalize the approach to the replicated 
model, EqO] allowing us to study the role of quenched dis- 
order. There is a close similarity of this approach to the renor- 
malization group developed for many body fermion systems, 
where the low energy modes are also located on a sphere of 
finite momentum, the Fermi surface i4lll . The distinct aiTange- 
ments of the scattering momenta lead in case of fermions to 
BCS, forward and exchange scattering. This similarity sup- 
ports that the Brazovskii model will be characterized not only 
by one four point vertex, but by an effective interaction with 
generic angle between incoming and scattered momenta and 
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an interaction where all interactions are parallel. 

The most general form of the quartic interaction in Eall His 
expressed in terms of two vertex functions F and A: 



qi,q2.q3.q4Vo,qi 



I ^(ii.,(i-,.,ii^.,(i3,^u.(\i^a,(\2^b.(\^^h,Hi, 



with shorthand T, 
the inter-replica coupling vertices are given by 



The initial values of 



Aq,,q2,q3,q4 = -AqS (qi+q2+q3+q4) . 

The vertex Tq^ q^ q^ q^ has been determined by Leiblerfl. 
Considering |q,| = it only depends on two independent 
angles 9i2, between vectors qi and q2, and 9i4 between vec- 
tors qi and q4, respectively. In the notation of Ref.|2], we 
can write rq(^q2,q3,q4 =T4{hi,h2) , where /ii =2 + 2cos(9i2) 
and h2 = 2 + 2cos (9 14). The most relevant vertices are those 
where all momenta lie in the same plane. This only happens 
if either hi —0 or /12 = (equivalent to 9i2(4) — Tl and 9i4(2) 
generic) or /ii +/12 = 4 (equivalent to 9i2 = 9i4 ±7i). We find 
that 



r4(cos(9)) = -^fl,P/(cos(9)), 



(12) 



where f, are Legendre polynomials and the coefficients of 
the first three angular momentum channel are oq — 176.92, 
fli = 14.34, 02 — 2.30 and 133 — 8.33. Furthermore, the angle 
9 stands for 9 = 9i2 or 9i4. In the course of the renormaliza- 
tion, these vertices may change differently depending on the 
specific momentum dependence. The vertex F, which does 
not couple distinct replicas, leads at low energies to the two 
coupling constants 



To = rp,_ 
roil = rp,_ 



p,q, 
p>p. 



where Fq refers to a generic angle between the momenta p and 
q (9i2 = ^ and generic 9i4) with |p| = |q| = qo, while Fo|| de- 
termines the renormalized interaction of the special case p = q 
(9i2 — 9i4 — k). As first shown in Ref.|8J, Fo and Fo|| renor- 
malize differently. Using EqElone can determine the initial 
values of these two coupling constants: 

106.20 

Toil = T^^, (13) 



Fo 



120.00 



(14) 



Naoc"^ ' 

The situation becomes more subtle in case of the vertex A 
which couples distinct replicas. The symmetry of this interac- 
tion allows for altogether four distinct vertices: 



A 


Ap, p.q, 


q 


A|| 


= Ap, p.p, 


p 


u 


= Ap q, p 


q 


^11 


= App, p, 


p 



The coupling constants A and Q. refer to generic angles be- 
tween the two scattered momenta, whereas they are paral- 
lel in case of A|| and f2||, respectively. While F and F|| 
characterize the interactions of the system, the coupling con- 
stants A, A|| , Q. and £1|| determine the evolution of the disorder 
distribution function. A and A|| corresponds to local disor- 
der, as characterized by the bare distribution function, Eq|8l 
f(x)f(x') oc 5(x — x'). On the other hand, Q. and stand for 
new non-local randomness generated as one considers the ef- 
fective low energy properties of the system. This last aspect 
becomes most evident by undoing the replica trick. Introduc- 
ing 

Gp q — Ap_q, p _q , 

the corresponding part of the replicated action is 

SA = jTy\ [ G(x + x',y + y')(j)„,x(|)«x4i.y'l>iy ' 

with real space Fourier transform G(x,x') = 
r^f^^i(px+qx')(3: ^^jj considered as the 

■' (27C) ^'^ 

result of a replica calculation where the initial disordered 
Hamiltonian was 



The nonlocal random function T (x, y) has a distribution func- 
tion 



1 



P[T] - exp (^-^y d'^xd'^x'd'^yd'^y' 

r(x,y)G-'(x-y,x'-y')r(x',y')) , 

with correlation function G(x — y,x' — y') = 
(x,y) r (x',y'). Thus, in case Q. or fi|| flow to large 
values, the system acquires a non-local randomness with 
typical length scale |x — y|,|x'— y'| ~ In this sense we 
obtain a renormalization of the distribution function of the 
problem. The initial distribution 

Go (x - y , x' - y') = 2Ao5 (x - y) 5 (x' - y') 

becomes broadened such that 

5(|x-y|)^/(^o|x-y|). 

and the disorder correlations become nonlocal. 



The clean case: Brazovskii transition 

Before we analyze the disordered case with altogether six 
distinct coupling constants we summarize the renormalization 
group approach to the clean case, Aq = 0. The derivation of 
the corresponding flow equations was given in Ref.|28] and 
we repeat the main steps in appendix A for completeness and 



5 



in order make the description of the repHcated theory more 
transparent. In addition we determine the actual phase bound- 
ary from the flow equations and compare the resuh with Bra- 
zovskii's original calculation|46|. As shown in appendix A, 
the flow equations of the clean model and after rescaling of 

2 

the coupling constants and mass term according to 7= -4tj^, 

71-A5 

rrii and r = are: 



It- An 



dr 
dl 



= 2r + 3y 



-dl - '^11^'^- 
This set of equations has a closed solution: 



r(/) 
Y(0 

Y|| (0 

where 7(0) ^% i 

_ to 



3 / dxyix), 
Jo 



-2x 



(15) 



(16) 



l+Yoil ie"-l) 



-yo\\e 



% 



ro,Yo|| =Y|| (0); 



Fo and r (0) 



ro = -4. Note that y(/) is always positive and approaches a 

fixed value y(/ qq) 1. However, Y|| (/) changes sign for 
I — I*, where /*is given by 



/* 



1 



■log 1 



1 

Yo 



(17) 



(more precisely at the point where metastable ordered states 
emerge). From Eqs^]and[n]the condition Iq — I* is obeyed 
if 



ro = e-^'* -3 dy 
Jo 



1 



We assume for simplicity Yo = Yo|| = Yo- In the limit Yo ^ 1 
holds e' — * Yo ''^^ ^"d we obtain to leading order in Yo 



V3 
1.507 Yf 



In (2) 



2/3 

Yo' - 



0{yo) 



(18) 



Returning to the unsealed variables and using the definition 
EqlSjfor To.we obtain that the stability limit (i.e. the spinodal) 
of the lamellar phase is given by 



(XA^), -10.50 - 



a 



(19) 



with a = 35.69. The value for a is comparable to the result 
a = 41, obtained in references H fioll using the self consis- 
tent Hartree theory of Ref.|8]. If we further assume the initial 
vertices according to EqO instead of Yo — Yo||' we obtain 
2/3 

— — 1.783 Yo , leading to a = 42.22. The remaining small 
difference is likely caused by the fact that our flow equations 
are only approximately valid for r(/o) ~ 1, no matter how 
small the coupling constant. Still the deviation in the numer- 
ical prefactor is very small and we will continue our analysis 
for the disordered case along similar lines. In addition we will, 
for the rest of the paper, assume Yo = Yo — Yo|| > with Fq given 
by EqEl 



If the system establishes a modulated order 

(x) = (|)o (x) exp {~iq,„ ■ x) , 

where |qo| = <?o and (|)o (x) is a smoothly varying function on 
the scale <7q the driving interaction is the one where all in- 
teracting momenta are either parallel or antiparallel, i.e. Y||- 
A flow towards negative Y|| indicated therefore a fluctuation 
induced first order transition to a modulated state, just like 
in Ref.[^. Since the scaling dimension of the interaction 
is three, we have, in distinction to a recent application to a 
quantum version of the Brazovskii modell'TT], no controlled 
e-expansion. Thus, we have to limit our analysis to the case 
of small Yo and Yo||- The RG-flow will proceed until scaling 
stops at a scale scale Iq with r (/q) = 1, i-e , To {h) — Aq. If 
Iq < I* scaling stops before the interaction changes sign. Then 
there is no first order transition. On the other hand, if Iq < I* 
the system changes character before scaling stops and we can- 
not determine the bare parameters without introducing a term 
of order (|)^. Thus we are in the ordered state or at least have 
local minima like in an overheated system. The transition hap- 
pens for Iq — I* and we are right at the first order transition 



The disordered case: flow equations 

Since all the six coupling constants listed above are allowed 
by symmetry, we have to start our calculation with a model 
where all those terms are included. We use the same rescaling 



of the coupling constants as in the clean case, i.e. 5 = 



Q., etc. The initial values of the flow are then 5 (0) = 



5|| (0) = o5(0) = (Oil (0) = -g, with 



Ao. 



(20) 



Generalizing the steps which led to the flow equations in the 
clean case to the replicated model Eq^Jleads to the following 
flow equations up to one loop: 



dr 



— = 2r + 3Y+2co 
dl 

^ = 3Y-(3f + 4Yra) 
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c/Y|| 

'dT 
S 

'dl 

'df 
dJ3 
'dl 

t/(0|| 

"df 





-2(3f 


+ 4703) 


35- 


- (6y5 + 
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4co5) 


3o|| 


— 1 oyo- 


-45© 4 


3co 






3cO| 







(21) 



Here the limit m ^ of the numbers of replicas was taken. 
This system of coupled flow equations can be solved in a 
closed fashion. The details of the solution are summarized 
in appendix B. The main results of this calculation are as fol- 
lows: For g < gc with 



3^/2(^/2- 



7(0 > Til (0 




1.4g. 
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■yo~0.22Yo, 



(22) 



0.5 1 1.5 2 / 



disorder does not change significantly the fluctuation induced 
first order transition, whereas for g > gc the system is domi- 
nated by strong, non local disorder fluctuations and no ordered 



this criterion can be ex- 



lamellar state forms. Since z- 

Yo 

pressed in terms of a critical value for the fluctuations of the 
Flory-Huggins parameter, y A%^, and is given by: 



(23) 



Disorder affects long polymer chains stronger than short 
chains, a consequence of the relation between the interaction 
and disorder strength given in EallOl 

If g < gc, the coupling constant y|| {I) changes sign (see fig- 
ure|3a) at a scale 



1, {,A^. 



(24) 



with 



1 - 2jt: + Vl 8x2 - 8x ■ 



No other coupling constant diverges or changes sign for I < 
I*. The behavior is similar to that of the clean Brazovskii 
transition and for ^ ^ we recover the clean limit, EgfTTI 
On the other hand, for g > gc all coupling constants diverge at 
the scale 



From the solution of Eq|2]it further follows that these diver- 
gencies are "driven" by the divergence of the coupling con- 
stant CO. As discussed in detail above, a large value of co 



Figure 2: RG flow of the coupling constants as derived from equation 
Ellfor Yo = 0.01 and g = 0.6gc, gc , l-4gc. a) Dashed (full) curve 
refers to y (y||). b) Dashed (full) curve refers to 5 (6||) 



implies a strong renormalization of the distribution function 
of the randomness. Nonlocal disorder correlation, which are 
clearly tied to strong random compositional fluctuations on 
the scale 2k /qo occur instead of an ordered lamellar state. 
In Fig. 2 we show the flow of the various coupling constants, 
demonstrating that for g < gc the coupling constant yi changes 
sign well before all other coupling constants diverge. On the 
other hand, for g > gc the systems flows to a behavior with 
strong randomness and nonlocal disorder correlations while 
y|| remains positive (see figure|2j-a). 

The details of the phase boundaries between the various 
states are determined by the flow behavior of the mass pa- 
rameter r{l). The solution for the mass flow equation is (see 
also figure|3j-a) 



ril) 



e^'[ro + e-^'' [2co(/') +3y(/')] 

^''('•o + S^yo-fjx 

e''dl' 



i + (yo-fj (.^''-1) 

Here, the divergencies of y and co at Ig are all canceled and the 
result for r (/) is precisely the same as for the clean system but 
with 

yo^yo-y. 

As in the clean system, scaling stops at r (/q) — 1 which leads 
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Figure 3: RG flows as derived from equation 12 II for yo = 0.01 and 
g = 0.6gc, gc ,i-^gc- a) for f-Q = rQ(see equation^}. Tlie dotted 
line defines tfie scaling limit Xq = Aq, i.e, r = 1. b) Dashed (full) 
curve refers to (3 ((B||) 
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to a boundary 



ro = e-'^'o - 3 



£•'0 



dy 



If ^ < gc — ^ g ^ Yo there is a fluctuation induced first order 
transition with finite randomness at Iq — I*, with /* defined by 
l24l On the other hand, for g > gc up to the leading order on 
g and Yo the boundary between the liquid and glassy phases is 
given by 



ro 



f) 



2/3 



271 / 

7=3^- 



3 



2/3 



In Figl^we show the phase boundary lines in the ro,Yo and g 
parameter space. 

A qualitative understanding of these results can be obtained 
by combining simple droplet arguments with a criterion simi- 
lar in spirit to the one developed by Harris f35*| for disordered 
second order phase transitions. We consider the equation of 
state of the lamellar phase in the clean system(3l, using the 
dimensionless variables of this paper: 



-'"0 



;-2 



Here, ^ is the dimensionless correlation length (the correlation 



length measured in units of '). This gives a value Tq = 

3^/3 

— where for the first time a nontrivial solution becomes 
possible (spinodal line). The value of the correlation length at 
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Figure 4: Phase boundaries of the model, Eq. |4|determined from the 
solution of the flow equations. Has. 1211 In the three graphics the solid 
lines represent the boundary between the liquid and lamellar phases, 
the dashed lines represent the boundary between the liquid and glassy 
phases and the dotted line the boundary between the lamellar and 
glassy phases. 

the transition is t,* — . In the presence of disorder, 

this can to leading order in Yo be generalized to 



-'"0 



-f(x) = r + Yo((l)x): 



where ((|)^) is the (appropriately rescaled) mean square devi- 
ation of (|)x in the presence of disorder t (x) is the random 
and has also been rescaled to be dimensionless. We consider 
a solution with fixed modulation direction qo (with |qo| — qo) 



along the JC|| -axis: (|)(x||,Xx) = A (xpx±^ e"^°^K As shown in 
Ref.|281, domain walls where qo || n, with domain wall nor- 
mal vector, n, are governed (3||A) and decay on the length 
scale of the correlation length, t,. On the other hand, do- 
main walls with qo -L n are determined by ^ (V^A)^ and 

decay on the scale yj^q^. For weak coupling ^qo ^ 1, and 
the transverse domain walls (qo _L n) of a droplet are much 
less costly than the longitudinal ones. Droplets will there- 
fore be ellipsoidal (cigar like). Disorder fluctuations along 
qo will predominantly affect droplet formation and are deter- 
mined by f|| (.\:||) = / c/^'^'xx? (x||Xx), independent on the dis- 
order variation perpendicular to qo- Its distribution function is 

The typical value of the disor- 



der in a regime of linear dimension L is then 



tL = - 

L 
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Figure 5: Phase boundaries between lamellar , glassy (nonlocal dis- 
order correlations) and disordered states for a disordered symmetric 
copolymer with A' =10"* monomers as function of the average Flory- 

Huggins parameter x/Xo ™d 'ts fluctuation, ^ = ^ ^ . The refer- 
ence units are respectively Xo > the critical Flory-Huggins parameter 
for the clean case (Aq = 0) as given bv ll9l and XjA' = 10.50. 



is then gives as f| = ^. Disorder 



The typical value Af^ 
is not changing the equation of state dramatically if At^ <^ 

which gives g ^ Due to the first order character of the 

— 1/3 

transition the correlation length, ^ never exceeds ^* ~ Yo ■ 
Thus, there is a regime at small disorder strength where the 
equation of state is not dramatically changed. The result- 
ing criterion is ^ < y in agreement with our renormalization 
group results. The limited correlation length at the first or- 
der transition protects the ordered state from disorder fluctua- 
tions. If the first order transition becomes too weak, the cor- 
relation length at the transition becomes arbitrarily large and 
the system becomes effectively critical. In this case the disor- 
der starts dominating the low energy physics and the system 
changes character. This conclusion is fuUy consistent with the 
renormalization group analysis presented above. It is the flow 
towards negative Y|| and the related first order transition, that 
avoids the otherwise inevitable disorder driven divergence of 
the coupling constants. 

To summarize our results we can plot a phase diagram (Fig- 
ure |5} in terms of the di-block copolymer parameters for a 
fixed number of monomers (A^ = lO'*) as a function of the rel- 
ative intensity of disorder. Ax = \J Ax^, as defined by equation 
|9] The boundary between the glassy state, defined bvl22land 
[T9lappears at ^ -0.1. 



CONCLUSION 

In summary, we have investigated the effect of a ran- 
dom Flory-Huggins parameter x (x) in a symmetric di-block 
copolymer melt. Such fluctuations occur in a copolymer 
where one block is near its structural glass transition! 26^ .27.1 . 
The decreased mobility of the monomers of this block im- 
plies a partial annealing and leads to slow non-thermal com- 
positional fluctuations L25il . In the clean limit the microphase 



segregation between the two blocks causes a weak, fluctua- 
tion induced first order transition to a lamellar state. Using 
a renormalization group approach combined with the replica 
trick to treat the quenched disorder, we showed that in case of 
small fluctuations of x (x) the first order transition to a lamel- 
lar state is unchanged by disorder. Once the strength of the 
spatial fluctuations of the Flory-Huggins parameter exceeds a 
critical value that depends on the length of the polymer chains, 
the character of the transition changes. The system becomes 
dominated by strong randomness. Very likely a glassy rather 
than ordered lamellar state occurs. In this disorder dominated 
regime a strong renormalization of the effective disorder cor- 
relations occur. Nonlocal disorder correlations emerge that 
reflect strong compositional fluctuation on the scale of the ra- 
dius of gyration of the polymer chains. If non-thermal sta- 
tistical fluctuation of the Flory-Huggins parameter within the 
volume occupied by a copolymer chain exceeds a fraction of 
the order of A^"'/'* of the critical coupling (w Xi — 10.49/A^) 
the lamellar state is unstable. Instead it leads to a state with no 
long range-order and a renormalized distribution of the Flory- 
Huggins constant. The strength of the first order transition 
depends on the chain length, and is more pronounced for 
shorter chains. In this case the correlation length of the system 
at the clean phase transition is comparatively short and criti- 
cal fluctuations are unimportant. The first order transition is 
then unaffected by weak disorder On the other hand, longer 
chains imply a much weaker first order transition and critical 
fluctuations with large characteristic correlation length come 
into play. These critical fluctuations are extremely susceptible 
with respect to disorder, which is reflected in a flow towards 
infinite disorder strength. The same result was obtained using 
a modified Harris criterionfss'l. The large phase space of the 
microphase separation transition in symmetric copolymers is 
therefore a very interesting example for the unique role that 
disorder can play at a first order phase transition. 
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APPENDIX A: THE CLEAN BRAZOVSKII MODEL 
WITHIN THE RENORMALIZATION GROUP APPROACH 

In this appendix we summarize the main steps of the renor- 
malization group approach to the clean Brazovskii model of 
Ref.j3- The bare action of this problem is 

_q. 



with q' = q + Qi + q2 and q" = q + qi — q2- In distinction 
to the usual (|)^-theory we need to consider these three differ- 
ent terms separately. There are two distinct ways to arrange 
the angle of the unit vectors ni and n2, leading to Y and r||. 
For a generic angle between q i and q2 only one of the above 
three terms has propagating lines which are always on the low 
energy surface, yielding 



-2 2 
Ja 



while in case of qi = q2 two of the three terms contribute and 
it follows 



^rqi,-q,,q,,-q, — 61^ / Gq. 



We introduce k = \a\ — and can write f 

111 J J ^2%f 

?o ' where dO. is the measure for the solid an- 

gle integration. 

A tree level analysis with k' — k/b and (|)^,, = b^^(^if gives 



1 

So = T^'^-' / dk'd£l(^',, ,, {k'^ + bh) c^'_,, ^ , 



Ao 
Ao 

yielding i' — b^X and ^ — \ for the scaling dimension of (|), 
independent of dimension d. 

In our treatment of the interaction term we also follow 
closely the approach of Ref. L41.I .The interaction term is given 

as 

~ A ^PlP2P3-(Pl+P2+P3) 
^ •'Pl,P2,P3 

Xe(A- |Pl +P2+P3|)(^pi(^p2(l)p3(^-(P1+P2+P3)' 

where the 0-function ensures that P4 = — (pi +P2 + P3) falls 
indeed into the small shell around the d — 1 dimensional sur- 
face of low energy excitations. In order to make the theory 
renormalizable the 9-function is softened via 9 (Aq — |p4|) 
exp(— |p4| /Ao). At the tree level it follows: 



-(6-l)<?o(|A|-l)/Ao;^3p 

" ^ k yb ■> 



with A= ni + n2 — n3, where n, are the unit vectors of the mo- 
menta p,-. The constraint |A| — 1 is a consequence of the fact 
that the momentum p4 lies within a shell [^0 ~ Ao,^o + Aq]. 
This implies that the four wave vectors p, are equal and op- 
posite in pairs when their magnitude equals ^o- Thus, only 
Fqi qj _qj survivcs and satisfies 



f: 



qi, qi,q2. q; 



2.-q2 (^qi -qi,q2,-q2 +^rq[,q2) ■ 



For a scalar field holds generally: 

Srqi,q2 = 



q>-q,qi -qi "q -q>q2 -q2 ^(fi<i 



3 r u 

Jq 

3 / Wq,-q',qi,q2Wq',-q,-qi,-q2^q^q' 
Jq 

3 / Un _n" n, -n^U„ll 
Jq 



qi -q2'*q " -q,q2, qi^q^q"- 



In both cases is the coupling constant in the interaction that 
for a generic angle, i.e. F. It holds 



Ja 



Ao 



dp 



1 



= -a- 
3 / 



AoZ 



.2N2 



where a : 
tions 



6qi-^Kd. Thus we find the following flow equa- 



dl 



3F|| 



gApF 



Similarly, it follows for the mass term t' = b^ {z + E), where 



3 / Tq/ _q/ q _q Gq, 
Jq 



which is only weakly momentum dependent and it holds 



dz 

— =2T + aAo— 
dl T- 



V2- 



Considering a weak first order transition we neglect 1 in the 
various denominators, i.e. t + Aq~Aq. We verified by nu- 
merically solving the full flow equations of the problem that 
this approximation only causes very minor changes in the fi- 



nal results. Rescaling y = and similarly for all other the 



3A; 



coupling constants via as well as r = yields the equations 
given in EallSI 



APPENDIX B: DETAILS OF THE SOLUTION OF THE 
FLOW EQUATIONS 

In this appendix we give the full flow equation for the disor- 
dered Brazovskii model and summarize in detail its solution. 
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A one loop calculation along the lines of the clean model, but 
taking into account all six coupling constants yields 



dl 
'dl 

~dr 

dA 
'dl 

dA\\ 

~dr 

dU 

'dl 
dQ. 

~di 



2'i + K{%) (3r + mA + 2n 
-2 



3r-n(T) (^3r +4raj 
3r|| -2n(T) (3r^+4ra) 

3A - n (t) [eTA + mA^ + AQZ}j 



^ 3A||- 



n (x) (erA + mt^ + \a£i + id 



3n - 2n (t) d 



—f- = 3n||-4n(T)i2 



(25) 



Since co (/ = 0) = —g < 0, co (Z) divergence at a value 

The important issue is whether the flow will ever reach Ig, i.e. 
whether there will be a divergence of other coupling constants, 
an instability of the system (Y|| < 0) or whether scaling stops 
for some / smaller than Ig. 

Considering next (M|| (/) leads to the solution (M|| (/) ~ 
2co(Z) + ge^' and v'|| (1) diverges together with v{l) and will 
not change sign for / < Ig- 

We next find the solution for y(/): 

10(5(1) /g 



W) = - 



1 



(yo-^)(^ 



I) 



which diverges if (m(/) does. If ^ > jjo, j{l) also diverges at 
the scale 



Here 



K{z) 

n(T) 



d 
dl 



1 Ao 
-a =■ 

3 T + A2 



1 



•■0 

Ao 



39? 



with a = We first approximate 11 (t) ~ n (0) and K (t) ~ 
K{0) which gives K{0) = and n(0) = We have 



solved numerically the flow equations without this simplifica- 
tions and only find very minor differences. Further we intro- 
duce dimensionless coupling constants Y = n(0)r and sim- 
ilarly for all other the coupling constants as well as r = -^j-. 

Performing the m ^ limit of the numbers of replicas finally 
yields the flow equations given in Eal^TI 

Before we give the solution of this set of equations we 
discuss the fixed points and demonstrate that no new dis- 
order fixed point occurs. There are altogether four differ- 
ent fixed points. At the Gaussian fixed point all coupling 
constants together with r* vanish. Next, we find the clean 
fixed point where all disorder coupling constants vanish, but 
Y* = 1,Y|=2;t* = — |. This fixed point corresponds to a 
tricritical point and was discussed in detail in Ref. ll47l . Fur- 
thermore, there are two more, unphysical fixed points where 
some of the co's and 5's are positive, or y's negative, which are 
in both cases signs that render the theory instable. Note, all 
the disorder variables have a bare value —g and are negative. 
Thus, there are no new disorder related fixed points. 

Next we summarize the solution of the flow equations. The 
above set of flow equations, Eq|21 can be solved in a closed 
fashion. It is useful to first solve the flow equation for (o(/); 



05(0 



•fie 



1 

lu = ^ log 



1 



However, for g > it always holds that Ig < l„, i.e. does 
not diverge before CO (l). 

Next we consider the solution for Y|| (0 which signals the 
fluctuation induced first order transition in the clean case. It 
holds 



Y|| il)=2y{l)-yoe" 

Y|| changes sign at the scale /* where j{l* 
above solution for y(0 it follows: 



70 



e . Using the 




l_il.±7 




with X {x) — V 1 + Sx^ — %x. As can be seen in Fig. 2, there are 
two /-values where y(0 changes sign. We have to look for the 
smaller value (the first sign change). The first time where this 
gives a real solution for I* happens when the argument of the 
square root in % [x] vanishes. This yields the g-value where 
Y|| changes sign first: 



3V2(V2 



- 1 



■Yo- 



Finally we discuss the solutions for 5 (/) and 5|| {I). It holds 



5(0 



(m_\ 



The solution, 5(0 will, for the same reason as y(0' ^lot di- 
verge before ro(/). Considering 5|| (0: 



„3Z 



§11 (0 = 



(l + ^(.^^-lf (^-Yo)y(O) 
(l-f (.3'- 1)) (1+3,(0) 
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with y{l) = (e^' — l) — Yoj — 2. 5|| diverges at the scale 
Ig. However for ^ > ^Jq the coupling constant a)|| (l) changes 
sign. Restricting ourselves to g < ^yo we do not need to 
specify whether this sign change ever takes place before c5 (/) 
diverges. 



[1] E. Helfand, J. Chem. Phys. 62, 999 (1975). 
[2] L. Leibler, Macromolecules 13, 1602 (1980). 
[3] T. Ohta, and K. Kawasaki, Macromolecules 19, 2621 (1986). 
[4] G. H. Fredrickson, E. Helfand, J. Chem. Phys. 87, 697 (1987). 
[5] F. S. Bates and G. H. Fredrickson, Annu. Rev. Phys. Chem. 41, 
525 (1990). 

[6] R. Holyst and T. A. Vilgis, Macromolecular Theory and Simu- 
lation, 5, 573 (1996). 
[7] 1. W. Hamley,r/je Physics of Block Copolymers (Oxford Uni- 
versity Press, Oxford) 1998. 
[8] S. A. Brazovskii, Soviet Phys. JETP 85 (1975). 
[9] S. A. Brazovskii and S. G. Dmitriev, Zh. Eksp. Teor. Fiz. 69, 
179 (1975); S. A. Brazovskii and V. M. Filev, Zh. Eksp. Teor. 
Fiz. 75, 1140 (1978). 
[10] F. S. Bates, J. H. Rosedale and G. H. Fredricksen, J. Chem. 

Phys. 92, 6255 (1990). 
[11] T. Garel and H. Orland, Europhys. Lett. 6, 597 (1988). 
[12] E. 1. Shakhnovich and A. M. Gutin, J. Phys. (Paris) 50, 1843 
(1989). 

[13] G. H. Fredrickson and S. T. Milner, Phys. Rev Lett. 67, 835 
(1991). 

[14] G. H. Fredrickson, S. T. Milner, and L. Leibler, Macro- 
molecules 25, 6341 (1992). 
[15] S. R Panukov and S. I. Kuchanov, Sov Phys. JETP 72, 476 

(1993). 

[16] 1. Ya. Erukhiomovich and A. V. Dobrynin, Sov. Phys. -JETP 

Lett. 57, 125 (1993). 
[17] A. V. Dobrynin and 1. Ya. Erukhiomovich, Sov Phys.-JETP 77 

307 (1993). 

[18] A. M. Gutin, S. D. Sfatos, and E. I. Shakhnovich, J. Phys. A: 

Math. Gen 27, 7957 (1994). 
[19] J. Houdayer and M. MuUer, Macromolecules, 37, 4283 (2004). 



[20] R G. de Gennes, J. Rhys. Chem. 88, 6469 (1984). 

[21] S. Stepanov, A. V. Dobrynin, T. A. Vigilis and K. Binder , J. 

Rhys. 1 France 6, 837 (1996). 
[22] J. Schmalian and RG. Wolynes, Rhys. Rev Lett. 85, 836 (2000). 
[23] H. Westfahl, Jr., J. Schmalian, and RG. Wolynes, Phys. Rev B 

64, 174203 (2001). 
[24] S. Wu, H. Westfahl, Jr., J. Schmalian, and RG. Wolynes, Chem. 

Phys. Lett. 359, 1 (2002). 
[25] A. V. Dobrynin, J. Phys. I France 5, 657 (1995). 
[26] M. V. Giotto, C. L. Sangiorge, D. J. Harris, A. L. deOliveira, 

K. Schmidt-Rohr, T. J. Bonagamba, Macromolecules 35, 3576 

(2002). 

[27] L. Zhu, Y. Chen, A. Zhang, B. H. Calhoun, M. Chun, R. R 
Quirk, S. Z. D. Cheng, B. S. Hsiao, F. Yeh, and T. Hashimoto, 
Phys. Rev B 60, 10022 (1999). 

[28] R C. Hohenberg and J. B. Swift, Phys. Rev E 52, 1828 (1995). 

[29] R Bak, S. Krinsky, and D. Mukamel, Phys. Rev Lett. 36 , 52 
(1976). 

[30] Y. Imry and M. Wortis, Phys. Rev B 19, 358 (1979). 
[31] M. Aizenman and J. Wehr, Rhys. Rev Lett. 62, 2503 (1989). 
[32] K. Hui and A. N. Berker, Rhys. Rev Lett. 62, 2507 (1989); ibid 
63, 2433 (1989). 

[33] J. Cardy and J. L. Jacobsen, Phys. Rev Lett. 79, 4063 (1997). 
[34] J. Cardy, Rhysica A 263, 215 (1999). 
[35] A. B. Harris, J. Rhys C 7, 1671 (1974). 
[36] T. C. Lubensky, Rhys. Rev. B 11, 3573 (1975). 
[37] D. E. Kmel'nitskii, Sov Rhys. JETP 41, 981 (1976). 
[38] G. Grinstein and A. Luther, Rhys. Rev B 13, 1329 (1976). 
[39] A. Aharony, Y. Imry, and S.-K. Ma, Phys. Rev B 13, 466 
(1976). 

[40] K. G. Wilson, Phys. Rev B 4, 3174 (1971); ibid 3184 (1971). 
[41] R. Shankar, Rev Mod. Phys. 66, 129 (1994) 
[42] A. V. Dobrynin, and Erukhimovich, 1. Ya., J. Rhys. 11 France, 1, 
1387 (1991) 

[43] Dobrynin, Audrey V. and Erukhimovich, Igor Ya.,J. Phys. I 

France, 5 , 365 (1995) 
[44] Dobrynin, A. V. and Stepanow, S. and Vilgis, T.A. and Binder, 

K., J. Rhys. 1 France, 6 , 837 (1996) 
[45] A. Weinrib, and B. I. Halperin , Rhys. Rev B 27, 413 (1983). 
[46] K. Bates, Adv Pol. Sci 112, 181 (1994) 
[47] J. Schmahan and M. Turlakov Rhys. Rev Lett. 93, 036405 

(2004). 



